The superconductivity of Sr2RuO4 under c-axis uniaxial stress

Applying in-plane uniaxial pressure to strongly correlated low-dimensional systems has been shown to tune the electronic structure dramatically. For example, the unconventional superconductor Sr2RuO4 can be tuned through a single Van Hove point, resulting in strong enhancement of both Tc and Hc2. Out-of-plane (c axis) uniaxial pressure is expected to tune the quasi-two-dimensional structure even more strongly, by pushing it towards two Van Hove points simultaneously. Here, we achieve a record uniaxial stress of 3.2 GPa along the c axis of Sr2RuO4. Hc2 increases, as expected for increasing density of states, but unexpectedly Tc falls. As a first attempt to explain this result, we present three-dimensional calculations in the weak interaction limit. We find that within the weak-coupling framework there is no single order parameter that can account for the contrasting effects of in-plane versus c-axis uniaxial stress, which makes this new result a strong constraint on theories of the superconductivity of Sr2RuO4.

c-axis upper critical field H c2 increases by a factor of twenty 29 . This very strong enhancement is qualitatively consistent with, for example, a d x 2 Ày 2 order parameter. Under compression along the [100] direction, the Lifshitz transition occurs at approximately (k x , k y ) = (0, ± π/a), which we label the Y point. As the transition occurs, the Fermi velocity at the Y point falls to nearly zero, which in general is expected to increase T c and H c2 of order parameters such as d x 2 Ày 2 in which the gap is large at the Y point. The data under [100] uniaxial stress argue against, for example, a d xy order parameter.
Naively, then, T c and H c2 might be expected to rise even further under compression along the c axis. c-axis compression raises the energy of the d xz and d yz bands relative to the d xy band, and the resulting transfer of carriers expands the γ sheet, pushing it towards a Lifshitz transition from an electron-like to a hole-like geometry 30 . This transition occurs at both the X and Y pointssee Fig. 1c so the increase in the Fermi-level density of states (DOS) as it is approached is expected to be larger than for the electron-to-open Lifshitz transition induced by in-plane stress. Under a-axis compression T c increases strongly well before the Lifshitz transition is reached, and so generically we expect this to occur for c-axis stress, too. The weak-coupling renormalization group study of Ref. 31 and functional renormalization group study of Ref. 32 both predict a rapid increase in T c with approach to the electron-to-hole Lifshitz transition.
The electron-to-hole transition has been approached, and crossed, in thin films through epitaxial strain, and in bulk crystals by substitution of La for Sr [33][34][35] , but in both cases the superconductivity was suppressed by disorder. Here, we apply up to 3.2 GPa along the c axis of Sr 2 RuO 4 . This is a record uniaxial stress for bulk Sr 2 RuO 4 , and was achieved by sculpting samples with a focused ion beam to concentrate stress. H c2 increases, as expected from the increasing density of states. However, unexpectedly, T c decreases. In other words, approaching the Lifshitz transition at either the X or Y point dramatically enhances T c , while approaching both suppresses T c . This is a major surprise. In a first attempt to address this issue we present calculations in the limit of weak coupling, that take into account the three-dimensional structure of the Fermi surfaces. Although these show that c-axis compression reduces the transition temperatures of certain order parameters, no order parameter could be identified for which the effects of both outof-plane and in-plane pressure were captured. Our experimental finding therefore consitutes a major new constraint on theories of the superconductivity of Sr 2 RuO 4 .

Electronic structure calculations
We start with density functional theory (DFT) calculations of Sr 2 RuO 4 under c-axis compression, as a guide to the likely effects of c-axis strain on the electronic structure. Figure 1  (To simplify discussion, we refer throughout this paper to the X and Y points defined in the 2D zone, rather than the points of the full 3D zone observed in ARPES measurements 36 . Technical details of the calculation are provided in the Methods section. The calculations predict that the electron-to-hole transition will occur at ε zz = −0.025. Under a-axis compression, these calculations predict that the electron-to-open transition occurs at ε xx = −0.0075, whereas it was observed experimentally to occur at ε xx = −0.0044 28 , so this c-axis prediction might similarly overestimate the level of compression required. The uncertainty arises from the fact that the distance to the Lifshitz transition is sensitive to meV-level energy shifts, likely driven by many-body renormalisation 28 . Low-temperature ultrasound data give a c-axis Young's modulus of 219 GPa 37 , so ε zz = −0.025 corresponds to σ zz ≈ −5.5 GPa. Separately, we note also that while k z warping increases on all the Fermi sheets, as expected for c-axis compression, the β sheet has the strongest k z warping both at ε zz = 0 and at the Lifshitz transition; see the Methods section for an illustration.

Experimental results
Four samples were measured. For good stress homogeneity, samples should be elongated along the stress axis, which is a challenge for the c axis because the cleave plane of Sr 2 RuO 4 is the ab plane. A plasma focused ion beam, in which material is milled using a beam of Xe ions, was therefore used to shape the samples. Sample 1 was prepared with a uniform cross section, and a large enough stress, σ zz = −0.84 GPa, was achieved to observe a clear change in T c . To go further, the other samples were all sculpted into dumbell shapes, with the wide ends providing large surfaces for coupling force into the sample. FIB microstructuring has been used to achieve large c-axis stress in CaFe 2 As 2 38 , but here we needed to retain sufficient sample volume for high-precision magnetic susceptibility measurements. For measurement of T c in the neck portion, two concentric coils of a few turns each were wound around the neck. Samples 1 and 4 also had electrical contacts, for measurement of the c-axis resistivity ρ zz . Photographs of samples 2 and 4 are shown in the Methods section.
Sample 4 was measured in apparatus that incorporated a sensor of the force applied to the sample 39 , from which the stress in the sample could be accurately determined. Samples 1-3 were mounted into apparatus that had a sensor only of the displacement applied to the sample, which is an imperfect measure of the sample strain because the measured displacement includes deformation in the epoxy that holds the sample. Therefore, a displacement-to-stress conversion was applied to samples 1-3 to bring the rate of change of T c over the stress range 0.92 < σ zz < − 0.20 GPa into agreement with that of sample 4. In other words, we impose on our data an assumption that the initial rate of decrease in T c is the same in all the samples, which is reasonable because their zero-stress T c 's are very similar: all are between 1.45 and 1.50 K.
We begin by showing resistivity data, in Fig. 2. The plotted resistivities are corrected for the expected stress-induced change in sample geometry (reduced length and increased width), using the lowtemperature elastic moduli reported in Ref. 37, and making the assumption that stress and strain are linear over this entire range. At zero stress the resistivity of sample 4 shows a sharp transition into the superconducting state at 1.55 K. This sharpness, and the fact that it only slightly exceeds the transition temperature seen in susceptibility, indicate high sample quality. With compression, T c decreases. The normal-state resistivity also decreases, following the general expectations that c-axis compression should increase k z dispersion.
We find elastoresistivities (1/ρ zz )dρ zz /dε zz , obtained from linear fits over the range − 0.5 < σ zz < 0 GPa, of 37 and 32 for samples 1 and 4, respectively. Sample 4 was compressed to −1.7 GPa, and its resistivity does not show any major deviation from linearity over this range. The scatter in the data at strong compression may be a consequence of cracking in the electrical contacts-we show below that the sample deformation was almost certainly elastic.
We now show the effects of c-axis compression on magnetic susceptibility. Figure 3a-c shows the transitions of samples 2-4 in susceptibility; the data shown are the mutual inductance M of the sense coils versus temperature. To check that sample deformation remained elastic, we repeatedly cycled the stress to confirm that the form of the M(T) curves remained unchanged; see the Methods section for examples. For samples 3 and 4, the transition remained narrow as stress was applied, indicating high stress homogeneity. For sample 2, there was a tail on the high-temperature side of the transition, that was stronger at higher compressions. We attribute it to in-plane strain, possibly originating in the fact that sample 2 was not as well aligned as samples 3 and 4. A similar, though weaker, tail is also visible for sample 3. We note that the width of the transitions in Fig. 3a-c -≈ 50 mKwill be a consequence of defects and/or an internal distribution in the in-plane strain. Although there will also be inhomogeneity in ε zz , this is not the driver of the transition width: the distribution would have to have a width of~1 GPa, which is not plausible. Figure 3 (d) shows T c versus stress for all the samples. T c is taken as the temperature where M crosses a threshold. For samples 1, 3, and 4, we select a threshold at ≈ 50% of the height of the transition, and for sample 2, 20%, in order to minimize the influence from the high-temperature tail. T c is seen to decrease almost linearly out to σ zz ≈ − 1.8 GPa. For sample 4 (to which, as described above, the other samples are referenced), dT c /dσ zz in the limit σ zz → 0 is 76 ± 5 mK/ GPa. The error is 6%: we estimate a 5% error on the calibration of the force sensor of the cell, and a 3% error on the cross-sectional area of the sample (155 × 106 μm 2 ).
At σ zz ≲ − 1.8 GPa, the stress dependence of T c flattens markedly. In sample 3, T c resumes its decrease for σ zz < − 3 GPa. We show in the Methods section that both the flattening and this further decrease reproduce when the stress is cycled, which, in combination with the narrowness of the transitions, shows that this behavior is intrinsic, not an artefact of any drift or non-elastic deformation in the system. Figure 4 shows measurements of the c-axis upper critical field. M(H) for samples 2 and 3 at constant temperature T ≈ 0.3 K is shown in panels a and b. In Fig. 4c, we plot H c2 versus stress, taking H c2 as the fields at which M crosses the thresholds indicated in panels a-b. H c2 increases as stress is applied, as generally expected when the density of states increases. The increase is faster for sample 2 than sample 3, which may be an artefact of the tail on the transition for sample 2.
The quantity H c2 =T 2 c is particularly informative: if pairing strength were modified without changing the gap structure, H c2 would be  proportional to T 2 c . As shown in Fig. 4d, we observe H c2 =T 2 c to increase by 40% by σ zz = − 3.0 GPa. Slower quasiparticles are less strongly affected by magnetic field, so what this result means is that the Fermi velocity decreases on portions of the Fermi surface where the gap is large.
We conclude this section with a note on the peak effectthe small maximum in the susceptibility just below H c2 , visible in Fig. 4 (a, b). The peak effect occurs when there is a range of temperature below T c where vortex motion is uncorrelated, allowing individual vortices to find deeper pinning sites 40 . The peak is suppressed by c-axis compression, and it is suppressed downward rather than by being smeared horizontally along the H axis, which means that this suppression is not an artefact of a spread in H c2 due to strain inhomogeneity. It could indicate stronger pinning, due to the reduction in the coherence length.

Weak-coupling calculations
As explained in the introduction, there is an apparent contradiction between the increase of T c under a-axis strain reported in previous work (which suggests anti-nodes at the X and Y points) and the decrease of T c under c-axis pressure (which suggests nodes at the X and Y points). To see if this puzzle has a straightforward solution, we perform weak-coupling calculations for repulsive Hubbard models, as developed in Refs. 41-50. To capture possible changes in the 3D gap structure, we employ three-dimensional Fermi surfaces 51 . These are described by a three-band (4d xy, xz, and yz) tight-binding model. The hopping integrals are derived from the Ru-centred Wannier functions obtained in the DFT calculation presented above. Our tight-binding model takes the form ψ s ðkÞ = ½c xz,s ðkÞ,c yz,s ðkÞ,c xy, s ðkÞ T , and H s ðkÞ incorporates spin-orbit coupling, inter-orbital and intra-orbital terms. The complete set of tight-binding parameters retained here is given in the Methods section.
In Fig. 5a, we show the tight-binding Fermi surfaces at ε zz = 0 and − 0.02. In Fig. 5b, we show the orbital weight on the γ sheet at k z = 0. As the γ sheet expands, the orbital mixing around its avoided crossings with the β sheet is reduced, and it becomes more dominated by xy orbital weight.
To H 0 we add on-site Coulomb terms projected onto the t 2g orbitals 52 (Methods Eq. (8)) and study the solutions to the linearized gap equation in the weak-coupling limit U/t ≪ 1, where U is the intraorbital Coulomb repulsion and t is the leading tight-binding term. We take the interorbital on-site Coulomb repulsion to be U 0 = U À 2J, where J is the Hund's coupling, and the pair-hopping Hund's interaction J 0 to be equal to the spin-exchange Hund's interaction J. Under these assumptions, the remaining free parameter is J/U. We take where μ and ν are band indices, |S ν | is the area of Fermi surface sheet ν, and Γ is the two-particle interaction vertex calculated consistently to order OðU 2 =t 2 Þ. Solutions to Eq. (2) with λ < 0 signal the onset of superconductivity, at the critical temperature T c~W expðÀ1=|λ|Þ, where W is the bandwidth.
In a pseudo-spin basis each eigenvector φ belongs to one of the ten irreducible representations of the crystal point group D 4h 48,55 . We calculate the leading eigenvalues in four even-parity channels, B 1g , B 2g , A 1g , and A 2gsee the legend of Fig. 5c-d. The E g channeld xz ± id yzhas been found to be strongly disfavored in weak-coupling calculations 51 , and so is not considered here.
A subset of the present authors have found, in a previous weakcoupling calculation, that the odd-parity order parameters track each other closely as J/U is varied, with a splitting that is small compared with that between the even-parity orders, and between the odd-and even-parity orders 51 . Ref. 56, likewise, finds the splitting between the odd-parity orders to be small. For this reason, we calculate only one odd-parity channel, E u (p x ± ip y ), and its behaviour can safely be taken to represent the qualitative behaviour of odd-parity order.
The leading eigenvalues in each channel as a function of ε zz are shown in Fig. 5c. Although, as in Ref. 51, odd-parity order is found to be favored, calculations in the random phase approximation at similar J/U tend to favor even-parity order 15,56 . A tendency towards odd-parity order appears to be a feature of calculations in the weak-coupling limit. We note also that the ordering of the channels differs from what was found in Ref. 51, due to a different tight-binding parametrisation. The ordering is sensitive to the parametrisation, and so we focus discussion here on trends with applied strain.
The weak-coupling results show a dichotomy in the strain dependence of T c : T c in the channels that have symmetry-imposed nodes at the X and Y points (E u , A 2g , and B 2g ) decreases with initial caxis compression. These nodes coincide with the regions of highest local density of states, and this result is an indication that order parameters in these channels are less able to take advantage of the increase in Fermi-level density of states induced by c-axis compression. However, under stronger compression T c increases modestly in all channels.
In the weak-coupling calculations of Ref. 29, the contrast in the response to a-axis uniaxial stress between order parameters with and without nodes at the X and Y points was found to be stronger in H c2 =T 2 c than T c , and so we also calculate H c2 =T 2 c , following the procedure in Ref. 29. Results are shown in Fig. 5d. We find that changes in H c2 =T 2 c correlate closely with shifts in the gap weight onto the γ sheet, which has the lowest Fermi velocity. For example, gap weight in the B 2g channel, shown in Fig. 5e, shifts from the β to the γ sheet as stress is initially applied, and H c2 =T 2 c correspondingly increases. At strong compression, gap weight shifts back to the β sheet, and H c2 =T 2 c decreases. This occurs because as the γ sheet expands it comes closer to its copies in adjacent zones, which disfavours a large gap on this sheet because in the B 2g channel the gap changes sign across the zone boundary. Among the even-parity channels, for ε zz < − 0.015H c2 =T 2 c increases for those without nodes along the Γ-X and Γ-Y lines, and decreases for those with. The complete set of calculated gap structures is shown in the Methods section.
We conclude this section by noting that although a k z dependence of the gap structure is seen in all channels, we do not find dramatic stress-induced changes in the k z dependence in any channel. Separately, in the A 2g channel there is a level crossing between ε zz = 0 and −0.0075. We plot only the leading eigenvalues in Fig. 5; this level crossing causes a large change in the leading gap structure and an anomalously large increase in H c2 =T 2 c .

Discussion
The unexpected decrease of T c as two Van Hove points in k-space are approached under c-axis compression is the key experimental result that we report. It might provide a vital clue about the nature of the superconducting state in Sr 2 RuO 4 , because it is so different to the response to in-plane, a-axis pressure. The DFT calculations indicate that our largest achieved stress, −3.2 GPa, is around 60% of the way to the Lifshitz transition, and if the calculations overestimate the   by its value at zero strain, of the leading eigenstate versus ε zz in each channel. In the legend, f i is any function that transforms as sin k i , d 0 (k) is the gap function for the even-parity irreducible representations, and d(k) the d-vector for E u . In the A 2g channel, there is a level crossing between ε zz = 0 and − 0.0075, that results an anomalously strong increase in H c2 =T 2 c . e Gap structure in the B 2g channel versus ε zz , on the α, β, and γ Fermi sheets. compression required to reach the Lifshitz transition, as they did for inplane stress, we might have come even closer. The contrast between the effects of a-and c-axis stress is unmistakable: 60% of the way to the electron-to-open Lifshitz transition under a-axis stress, T c is 0.7 K higher than in the unstressed sample 28 .
The response of T c under c-axis compression allows resolution of the stress dependence of T c into components through comparison with the the effect of hydrostatic compression, which also suppresses T c . We obtain the coefficients α and β in the expression where ΔV/V = ε xx + ε yy + ε zz is the fractional volume change of the unit cell, and ε zz − (ε xx + ε yy )/2 is a volume-preserving tetragonal distortion.
Refs. 12,57,58 report dT c /dσ hydro = 0.22 ± 0.02, 0.24 ± 0.02, and 0.21 ± 0.03 K/GPa; we take dT c /dσ hydro = 0.23 ± 0.01 K/GPa. Employing the low-temperature elastic moduli from Ref. 37 to convert stress to strain, we find α = 34.8 ± 1.6 K and β = − 2.2 ± 1.2 K. (Under hydrostatic stress, σ zz /ε zz = 396 GPa and ε xx = 0.814ε zz .) The small value of β means that a volume-preserving reduction in the lattice parameter ratio c/a would have little effect on T c : the increase in density of states by approaching the electron-to-hole Lifshitz transition is balanced, somehow, by weakening of the pairing interaction. The challenge for theory is to understand why that weakening takes place.
In the three-dimensional weak-coupling calculations presented here, it is the A 2g and B 2g channels, both of which have nodes along the Γ-X and Γ-Y lines, that best match observations. Due to differences between the actual and tight-binding electronic structure the ε zz = 0 point in the calculations should not be considered too literally as equivalent to ε zz = 0 in reality, and the key point is that it is only in the A 2g and B 2g channels that T c is found to decrease and H c2 =T 2 c to increase over some range of strain. However, as we have noted, A 2g and B 2g order parameters do not appear to be consistent with data under aaxis stress.
In other words, weak-coupling calculations do not explain the contrasting responses to a-versus c-axis stress, and this provides an opportunity: models of pairing in Sr 2 RuO 4 should be tested against this feature, for it might provide substantial resolving power between different models. There may, for example, be stress-driven changes in the interactions that drive superconductivity, though to attempt to calculate this is beyond the scope of this paper.
We highlight two other possible explanations. One is interorbital pairing [22][23][24][25] . The superconducting energy scale is too weak to induce substantial band mixing, and so these models depend on the proximity of the γ and β sheets, and the resulting mixing of xy and xz/yz orbital weight over substantial sections of Fermi surface 23 . We have noted that c-axis compression reduces this mixing, by pushing the γ and β sheets apart, which could then suppress T c 59 . In contrast, under in-plane uniaxial compression these sheets are pushed closer together along one direction and further apart along the other 36 .
The other is three-dimensional effects. Another feature of the electronic structure that varies oppositely under a-versus c-axis compression is the interlayer coupling. Under a-axis compression, the RuO 2 sheets are pushed further apart, and under c-axis compression, closer together. For example, an increase in warping of the Fermi surfaces along k z under c-axis compression could reduce the quality of nesting and so weaken spin fluctuations in Sr 2 RuO 4 , and the weakcoupling calculations here might not have fully captured the effect on the superconductivity.
In summary, we have demonstrated methods to apply uniaxial stress of multiple GPa along the interlayer axis of layered materials in samples large enough to permit high-precision magnetic susceptibility measurements. Under such a compression, we find that T c decreases even though the Fermi-level DOS increases, in striking contrast to the effect of in-plane uniaxial stress. Weak-coupling calculations do not provide a clear answer to this puzzle, which makes it important for models of superconductivity in Sr 2 RuO 4 to be tested against application of both types of stress. At a more general level, our findings motivate the use of out-of-plane stress as a powerful tool for investigation of other low dimensional strongly correlated system in which the strength of the interlayer coupling is suspected of playing an important role in their electronic properties.

Density functional theory calculation
DFT structure calculations were performed using the full-potential local orbital FLPO 60,61 version fplo 18.00-52 (http://www.fplo.de). For the exchange-correlation potential, the local density approximation applying the parametrizations of Perdew-Wang 62 was chosen. Spin-orbit coupling was treated non-perturbatively by solving the four-component Kohn-Sham-Dirac equation 63 . To obtain precise band structure and Fermi surface information in the presence of a Van Hove singularity close to the Fermi level, the final calculations were carried out on a well-converged mesh of 343,000 k points. (70 × 70 × 70; 23,022 points in the in the irreducible wedge of the Brillouin zone). As a starting point, for the unstrained lattice structure the structural parameters at 15 K from Ref. 64 were used. Longitudinal strain ε zz is taken as the independent variable, and ε xx and ε yy are set following the low-temperature Poisson's ratio from Ref. 37, which is 0.223 for stress along the c axis. The apical oxygen position was relaxed independently at each strain, by minimising the force to below 1 meV/Å. However, the effect of relaxing this internal parameter is small in comparison with the effect of the stress-driven change in lattice parameters.
The calculated Fermi surfaces of unstressed Sr 2 RuO 4 and under interlayer compression, including the warping along k z and the Fermi velocities, are shown in Fig. 6. The β sheet is the most strongly warped both at zero stress and at ε zz = −0.025.

Experimental details
Sr 2 RuO 4 samples were grown using a floating-zone method 65,66 . The four samples here were taken from the same original rod, and from a portion that we verified to have high T c and a low-density of Ru inclusions; our aim in taking multiple samples was to test reproducibility in sample preparation and mounting.
Uniaxial stress was applied using piezoelectric-driven apparatus 39,67 , and precision in sample mounting is important because Sr 2 RuO 4 is much more sensitive to in-plane than c-axis uniaxial stress: T c decreases by 0.13 K under a c-axis stress of σ zz = − 3.0 GPa, but increases by 0.13 K under an in-plane uniaxial stress of only 0.2 GPa 29 . Applying c-axis pressure could generate in-plane stress through bending and/or sample inhomogeneity. In a previous experiment 68 , c-axis compression raised T c and broadened the transition. However, the stress was applied at room temperature, where the elastic limit of Sr 2 RuO 4 is low 39 , so these effects may have been a consequence of in-plane strain due to defects introduced by the applied stress. Samples 2-4 were mounted into two-part sample carriers; that for samples 2 and 3 is diagrammed in Fig. 7a. The purpose was to protect samples from inadvertent application of tensile stress. Samples are mounted across a gap between a fixed and a moving portion of part B of the carrier, and can be compressed, but not tensioned, by bringing part A into contact with part B. In Fig. 7b, we show T c of samples 2 and 3 versus applied displacement, and the point where parts A and B come into contact and T c starts changing is clearly visible. For sample 2 the point of contact is rounded on the scale of a few microns, due to roughness and/or misalignment of the contact faces, and in all figures below we exclude data points that we estimate to be affected by this rounding.
The samples were mounted with Stycast 2850. This epoxy layer constitutes a conformal layer that ensures even application of stress 67 . Photographs of samples 2 and 4 are shown in Fig. 7c and d. The carrier for sample 4, which has a different design to those used for samples 2 and 3, is shown in Fig. 7e. Where electrical contacts were made, Du-Pont 6838 silver paste annealed at 450 ∘ for typically 30 min was used. This is longer than usual, in order to penetrate a thin insulating layer deposited during the ion beam milling.
As noted above, samples 1-3 were mounted in apparatus that had a sensor only of the displacement applied to the sample, while for sample 4 there was also a force sensor. Displacement sensors are less reliable as sensors of the state of the sample, because they also pick up deformation of the epoxy that holds the sample. In Fig. 7f the complete set of measurements of T c of sample 3, plotted against applied displacement, are shown. Data points are colored by the order in which they were collected. The data drifted leftward over time: stronger compression was needed to reach the same T c . However, the qualitative form of the curveinitial decrease in T c , then a flattening, and then further decreasereproduced over multiple stress cycles, and in Fig. 7g it is shown that the form of the transition was the same before and after application of the strongest compression. (We attribute the small apparent shift in T c to an artefact of inadvertent mechanical contact between the stress cell and inner vacuum can of the cryostat).
We therefore conclude that the sample deformed elastically and that it was the epoxy holding the sample that was compressed non-elastically; plastic deformation has previously been observed to broaden the superconducting transition 69 of Sr 2 RuO 4 . In Fig. 3, in the main text, we show only the data taken after the epoxy was maximally compressed. Force versus displacement data for sample 4 are shown in Fig. 7h-i, and here it can be seen that there was very substantial non-elastic compression of the epoxy. As with sample 3, the shape of the superconducting transition in the Sr 2 RuO 4 was the same before and after application of large stress. Over regions where the sample and epoxy deformed elastically, the combined spring constant was 1.45 N/μm. The spring constant of the flexures in the carrier, on the other hand, is calculated to be~0.03 N/μm, meaning that almost all of the applied force was transferred to the sample.

Calculated gap structure in other channels
In Fig. 8   A 2g channel, this shift occurs as a first-order change in gap structure between ε zz = 0 and ε zz = − 0.0075. At the largest compression reached, gap weight in the A 2g channel shifts back away from the γ sheet, as it does in the B 2g channel. This does not occur in the A 1g and B 1g channels.
For the interactions we use the (on-site) Hubbard-Kanamori Hamiltonian where i is site, a is orbital, and n ias = c y ias c ias is the density operator. We further assume that U 0 = U À 2J and J 0 = J 52 . In the weak-coupling limit this leaves J/U as a single parameter fully characterizing the interactions.
In the linearized gap equation (2) the (dimensionless) two-particle interaction vertex Γ is defined as 48 where ρ μ = |S μ |=½ v μ ð2πÞ 3 is the density of states, and 1= v μ = R S μ dk= |S μ |v μ ðkÞ . Here, Γ is the irreducible two-particle interaction vertex which to leading order retains the diagrams shown in Fig. 9. An eigenfunction φ of Eq. (2) corresponding to a negative eigenvalue λ yields the superconducting order parameter In the chosen pseudo-spin basis each eigenvector φ belongs to one of the ten irreducible representations of the crystal point group D 4h 48,55 .

Data availability
The data that support the findings of this study are openly available from the Max Planck Digital Library 70 .
performed the density functional theory and H.S.R., T.S. and S.H.S. the weak-coupling calculations; F.J. and C.W.H. analyzed data; and F.J. and C.W.H. wrote the paper with contributions from the other authors.

Funding
Open Access funding enabled and organized by Projekt DEAL.